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I. INTRODUCTION 

According to the paper by Nambu and Jona-Lasinio ^ or more precisely the Goldstone theorem, Q if the Lagrangian 
without scalar fields has a continuum global symmetry and if this symmetry is spontaneously broken, then there must 
appear dynamically generated massless spin-0 partilcles, the so-called Nambu-Goldstone (NG) bosons, the number 
' of which is determined how the symmetry is broken. Then it is necessary to have effective interations among the 
^ NG bosons and elementary particles to study the effects of this symmetry breakdown, either they appear as physical 
particles with tiny masses or are absorbed by guage fields to make them massive. 

Quite often these NG bosons acquire tiny masses through quantum corrections and/or soft symmetry breaking to 
Qv^ [ give actual physical effects. In these cases we really need an effective Lagrangian to calculate those physical effects. 

There are a couple of examples of effective Lagrangians known among the NG bosons and elementary particles. 
0^ ^ For instance, there are the original Nambu- Jona-Lasinio model Q of four-fermi interactions with spontaneous chiral 
0^ ■ symmetry breakdown in which the NG boson is completely massless, and the axion model |^ due to breakdown of the 
' Peccei-Quinn symmetry [Q in which the axion (pseudo-NG boson) acquires a little mass through quantum corrections, 
. etc. The former case, which is the first paper to stimulate the study of spontaneous symmetry breakdown, we just need 
O ^' to calculate the ordinary Feynman rules for elementary fermions as well as would-be NG bosons which are expressed 
D . as a chain of fermion loops. In the latter case we must assume and/or guess effective interactions with some unknown 
• coupling constants although some constraints are enforced. 

There is another interesting example in which people have calculated physical quantities without an effective action. 
When one calculates the pion decay constant, one uses the effective coupling between the pion (NG boson) and quarks 
(or nucleons), which was, utilizing the Jackiw- Johnson |^ sum rule, originally derived by Pagels and Stokar Q applying 
the Ward-Takahashi (WT) identity to the axial vector vertex. This way of determining the effective coupling, i.e., 
direct use of the WT identiy to the vertex, includes ambiguities. The Kyoto group has used the Bethe-Salpeter 
(BS) equation to uniquely determine this effective coupling and to calculate it perturbatively. They calculate 
the expectation value of an axial vector current sandwitched between vacuum and pion state and relates it with 
the BS amplitude to obtain the pion decay constant. In these approaches there do not appear the dynamical NG 
bosons. Instead Pagels and Stoker use only the WT identiy, and the Kyoto group uses a classical equation for the 
BS amplitude. Hence it does not seem to be clear how actually the NG boson couples to other elementary and/or 
dynamical fields as an operator form. 

Here we propose one approach to construct an effective Lagrangian among the NG bosons and elementary fields, 
i.e., the auxiliary field method, which was once in fashion to construct various kinds of models out of four-fermi 
interactions and also to study chiral symmetry breakdown. |^,^ With this method one can determine uniquely the 
coupling of the NG boson with fermions as an operator form. In this paper to show how powerful this approach is, 
we start from the simplest example, U{1) massless gauge theory without scalar fields to derive the coupling among 
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the dynamically generated NG bosons and fermions, to prove masslessness of the NG boson, and to calculate the 
decay constant of this NG boson in Section Then in Section Q] we proceed into SU{3) color gauge theory with 
iso-doublet fermions (up and down quarks) without scalar fields to do the same thing in Section || and to give the 
pion decay constant. The final Section ^ is devoted to the summary and discussions of our method compared with 
others. 



II. f/(l) MASSLESS GAUGE THEORY 



The Lagrangian for massless U{1) gauge theory is given by 
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(1) 



where = 9^ — igA^ and a is a gauge parameter. In this paper we use Minkowski metric except for a few equations. 
In order to introduce auxiliary fields as dynamical ones, we first functionally integrate the partition function over 
gauge fields A^. The resultant four-fcrmi interaction terms are rewritten by using the Fierz transformation as 

f d^xi^ix) %p^{x) + \g^ f d^xd^y^ix)^f,ijix)D'^''{x-y)ij{yhMy) 



d'^x ipix) ip^(x) 



d'^x d'^y 



'ij;{x)^{y)D{x - y)ip{y)'ijj{x) 



+ i^ix) ijb '>IAy)D{x - y)i^{x) 175 
where . . . expresses other modes and 



Df,^{x -y) = 
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(3b) 
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Then by introducing the gaussian term which is quadratic in auxiliary fields and cancels four-fermi terms in Eq. 



we obtain g, 
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d^x d'^y y>{x, y) - 'ip{x)ip{y)D{x - y)] D~'^{x- y) 
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X [7r{y, x) - D{y - x)i){y) 175 tpi^x)] + . . . , 



+ S'af = / d'^x ip{x) ipip{x) - I d^xd'^yip{x) <j>{x,y) + ij5Tr{x,y) ip{y) 



d'^x d^y 
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+ 
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(5) 



where D~^{x — y) is nothing but 1/D{x — y) and assumed are (l){x,y) = (f){y,x) and Tr{x,y) — TT{y,x). Since main 
purpose of this paper is to present our idea as simple as possible, we have neglected other modes. 

Next we assume that bilocal fields can be decomposed into local fields with continuum indices as follows, which is 
essential to the following whole arguments in this paper. 
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(X, y) = 0o(a; - y) H a 



iTQ{x-y) (x + y 



/ 



(6a) 
(6b) 



where fields with arguments x — y are considered to be classical fields and those with (x + y)/2 to be quantum ones like 
in the operator product expansion. Also assumed is that (j)o{x — y) is a vacuum expectation value (VEV) of (f){x,y). 
The VEV should depend only on x — y because of translational invariance. The coordinate {x + y)/2 is the center of 
mass system. Here a constant parameter / has massive dimension one, 4>o{x — y), <T{){x — y), and ^^{x — y) dimension 
five, and u{{x + t/)/2) and Lp{{x + t/)/2) dimension one. The difference between quantum and classical fields can be 
seen by Fourier-transforming fields with arguments x — y into momentum space, in which case, e.g., n^x^y) can be 
interpreted as local fields with continuum label of internal momentum. That is. 



TT{x,y) = 



T^o{x-y) fx + y 



f 



Tro{q)(p{X), 



with r — X — y and X ~ {x + y)l2, which shows 7r(a;, y) is equivalent to TTa{q) ip {X) with a continuum label q. 

The legitimacy of this decomposition Eqs. (^ may be supported by the following observation. The composite field 
in question is tightly bound due to the strong coupling which causes chiral symmetry breakdown and hence it can be 
decomposed into two wave functions as the first approximation, one for internal motion and another for total motion. 
An internal motion can be described in x — y, difference of coordinates of two particles, and the whole motion in 
{x + y)/2, the center of mass system coordinate of two particles. Anyway this is the simplest decomposition of (j){x,y) 
and 7t{x, y). This could also be justified by checking whether mass for the would-be NG boson, (p{(x + y)/2) becomes 
massless or not which we will see soon. 

From symmetrical point of view (rotational invariance) , it holds 



(^oix - y) TToix - y). 



(7) 



Since we are interested in describing a coupling of the NG boson, i.e., ^{{x -f y)/2), with fermions and in obtaining 
the equations for 4>q{x — y) and ttq{x — y), functionally integrating over fermions and keeping only 4>q{x — y), 7ro(a; — y), 
and ip{{x ~\- y)/2), we obtain 
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where A is given by Eq. ( ^ ) , Tr is to take trace over coordinates (a; and y or r — x — y and X = {x + y)/2 or X and q 
where g is a conjugate momentum to r) and gamma-matrix indices, is a dalambertian in terms of q^^ in Minkowski 
space and 



d^r 



Mr) e-'"^, Mq) 



d^r 



Mr) e"'"^ 



(9) 



In calculating Tr In, we have negelected the term inside since we are not interested in the derivative terms in X . 

Now setting ip{X) = in Eq. (H) since the VEV of this field is zero by definition {(p{X) is the NG boson) and 
varying ^off in terms of ^(9), we can obtain an equation for £((7). This is nothing but the gap equation, i.e., this 
gives the mass, ^(9), for a fermion. 
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or 
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(27r)2 (p-g)2p2 _5]2(p)' 



(11) 



which becomes a familiar form when converted into one in Euchdean space and in the Landau gauge a = or 
A = 3g^/{16iT^). Here we have used 



^ J (27r)2(p-, 



(27r)2(p-q)2- 

Next extracting only the bilinear terms in ip{X) from Eq. (P), we obtain the mass term for ip{X) as 

1 



2A/2 



d^x-^y^{x)fMq) 



92 + 4»A 
" ^ g2 _ s2(q) 



7ro(g), 



(12) 



(13) 



which is identically zero when one identifies T^oiq) with T,{q) and uses Eq. (p^. In other words, when one requires 
masslessness of the NG boson ip{X), then 7ro(g) = ^{q) is automatically derived up to a constant /. This proves 
masslessness of the NG boson, (p{X). Back to Eq. (||) this result means that the Yukawa coupling of the NG boson 
^p{X) with fermions is given by 



d xd yipix) i-f5 n (x, y) 



d p d q -~ i ( 1\ , \ ~ , \ 

— ;4 77r:^^(p + '?)j75S \v^T,) v{m\p). 



(14) 



in momentum space. Here fields with tilde are Fourier-transformed ones. The Feynman rule for the vertex of ipipip, 
i.e., fermion-anti-fermion-NG boson vertex, is finally given by 



This corrcponds to the Feynman diagram depicted in Fig.l with a vertex factor given by Eq.(p^. 



(15) 




Fig-l 



FIG. 1. Bare vertex of ij^ipxp. 



This result coincides with that derived by using the BS equation and slightly differs from that originally derived 
by Pagels and Stokar p] using the WT identity. Eqs. (o) are rewritten as 



4>o{x - y) 
f 



x + y 



<j)o{x~y) fx + y 
T^{x,v) = J V 



(16a) 
(16b) 



These equations mean that the fields (/)(x, y) and 7r(a;, y) should be rescaled by a factor {x — y) / f to obtain ordinary 
dynamical fields and a scalar field (j'{X) (= / 0(x, y)/(f>Q{r)) has a vacuum expectation value / to become a massive 
scalar field cr{X). 

Now we calculate the decay constant from this vertex. The decay constant in the U{1) case is given by 2/ and is 
defined by a matrix element of the axial vector current as 



(0|^(0) 7^75^(0)1 9^(9)) = -2z/9^. 



(17) 
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In the Qfi ^ limit, our vertex, the Beth-Salpeter vertex, and the Pagels-Stokar vertex, all give the same expression 
for the decay constant /. The factor 2 of Eqs. is explained in the final Section [V, wich is peculiar to U{1) 
theory. The Feynman diagram which corresponds to the left hand side (lhs)of Eq. (|17|) is depicted in Fig.2 which 
gives, keeping only a linear term in momentum and with the help of Eq. (|l5|), the final answer as. 



2 {27Ty 



xdx ■ 



Y.{x) Y.{x) ~xY.'{x)/2 



x + T?{x) 



(18) 



where momentum space is converted into Eucledean space, i.e., x = —p^. Be careful when one converts the results of 
Fig.2 to the r.h.s. of Eq. (p7| ) bacause is an outgoing momentum. 




FIG. 2. Fermion one loop correction to the axial vector current. 



III. SU{3) COLOR GAUGE THEORY WITH ISO-DOUBLET 

Now that we have described our idea in the former section, we would like to proceed into a more realistic case, 
SU{3) color gauge theory with massless iso-doublet fermions, u and d quarks. 

The difference between C/(l) and SU{3) cases is first the color factor, and secondly the flavor number. Although 
the arguments to obtain the pion decay constant in SU(3) are almost parallel to the U{1) case, there are a couple 
of numerical factor differences. Hence only the equations different from U{1) case are described below. First we have 
the following coupling beween the dynamical scalar and fermions. 



d'^xd'^yip{x) 0(a;,y) +i757r"(x,y)r° ^p{y), 



(19) 



where tpix) = (u(x), d{x)) and r" are the 2x2 Pauli matrices. Likewise we can easily derive masslessness of the NG 
boson, y)"(X), and the auxiliary fields given by Eqs. (^ should be replaced with 



TT°'{x,y) 



(j^ojx - y) 



x + y 



h{x-y) afx + y 



(20a) 
(20b) 



with 



d'^r 



:(j)Qir) e 



-iqr 



S(g). 



Instead of Eq. (fiTh for the U{\) axial vector current, the 5t/(3) axial vector iso-doublet current satisfies 



(21) 



^(0)7p75yV'(0) 



(22) 



Note that the coefficient of the rhs of Eq. (^2|) is 1 instead of 2. The Feynman rule for the vertex ij;ip°-ip in momentum 
space is given by 
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and we only give the final expression for as 



xdx ■ 



x + T?{x) 



(23) 



(24) 



where Nc = 3 is a number of colors. Only be careful that the Fierz transformation of internal degrees of freedom, i.e., 
SU{2) flavor as well as SU{3) color must be taken into account, in which case use has been made of 



(25a) 
(25b) 
(25c) 



<5a/3<557 = j;^Sa^Ssi3 + 2(r")^^ (T"')^/?' 

where T° are generators and iV is a number of dimensions of the group. 

IV. SUMMARY AND DISCUSSION 



In this paper we have shown how useful and easy our approach is and derived that by introducing the auxiliary 
fields, composite fields have a coupling with fermions and can be decomposed into local fields multiplied by a factor 
which is given by the solution to the gap equation, that the NG bosons are acually massless, and that the decay 
constant / is a vacuum expectation value of a composite scalar field and calculated by a one-loop fermion diagram 
given by Fig. 2 after taking a limit of 0. These have been done both for U{1) and SU{3) massless gauge theories 
without scalar fields. 

Now in this section, we show the relation in the simplest U{1) case between our method without the WT identity 
and others using the WT identity. If the chiral symmetry is not broken, the axial vector Ward-Takahashi identity is 
given by 



g^r^ (p + q,p)^Sp\p + q) 75 + l^Sp^ (p) , 
where Sf{p) is a fermion full propagator given by 

Sf{p) 



i> - E(p) ■ 

The usual way to obtain Eq. (p^) is to calculate the divergence of the T product. 



d^^e-^i-wv Q, (o| T {Jl{x)^p{v)m) |0) , 



J^ix) = i/'(a;)7A*75V'(a;) = - 



5C 



Sip. 



(26) 
(27) 

(28a) 
(28b) 



Here J^{x) is the Noether current for the chiral symmetry. Eq. (|28a| ) can be estimated by using the following 



anti-commutation relation, transformation laws, and conservation equation for the chiral symmetry, 



{ip\x,x°),my")} 



iS'^ix-y), 



[g^V(x)] =ij5iP{x), [Q^ipix)] =zV;(x)75, 



d^xJ^{x), 



d'^4{x)^Q, 



(29a) 

(29b) 
(29c) 



where is an operator which generates chiral transformation. Negative sign of Eq. ( |28b| ) comes from Eqs. ( |29I| ). 
Now let us graphically calculate T^{p + q,p) when the chiral symmetry is spontaneously broken. In this case, since 
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only the fermions as well as gauge fields are elementary, a massless pole contribution to r^(|j + q,p) can be depicted 
in Fig.3 and the whole T^^Xp + Q^p) with a tree vertex is given by. 



(p + p) = i 7m 75 - 2i / G{p + q,p)^ + regular part 
G(p + g,p) = is(p+|) 75, 
where "regular part" means the remaining terms which are not singular at = 0. 



(30a) 
(30b) 



7„7. <>C>C 




Fig.3 



FIG. 3. Massless pole correction to the axial-vector vertex. 



These are the argume nts o f References ^-0. Pagels and Stokar |^ have suggested the form of G{p + q,p) so that 
the WT identiy Eqs. (^6, 30a) is satisfied and gives a form slightly different from Eq. ( |301:| ) . The Kyoto group uses 
Eq. (30b) as a normalization of the BS amplitude. All these arguments have meaning o nly i n the limit of — > 0. We 
would like to put stress upon again the reason why the second term in the rhs of Eq. (30a) is added to F|^(p + q,p). 
This is because people so far have treated the NG boson as a composite particle of fermions, not an elementary field, 
and hence the Feynman diagram for the vertex shown in Fig.3 consisting only of fermions and gauge fields should be 
included. 

On the other hand, in our case we treat the NG boson as an elementary field although dynamical so that we can 
directly derive the coupling between fermions and the NG boson. Let us see how this affects the above arguments on 
derivation of the WT identity. Since the NG boson as well as a scalar field are introduced in the effective Lagrangian, 
the Noether current must be modified so that it includes the contributions from the scalars, which is given by 



sc 



5ip - 



SC 



(31) 



where J^(a;) is given by Eq. (28b). In our case the correct Noether current is given by Eq. (|3l|) and hence replacing 



Jfiix) with Jnix) in Eqs. (^91), the following equations hold. 



Q^V(a;) =^75^^(2;), Q^V(a;) 



a'ix) =2ip{x), Q^^ix) 



iip{x)-fr^, 
-2a'{x), 



d^xJ^ix), d^'J^Jx)=0. 



(32a) 
(32b) 
(32c) 



Using these equations, one can easily show that the same WT identity as Eq. (26) holds since Eqs. ( p2b| ) do not affect 
derivation of this equation. 

The second and third terms of the rhs of Eq. (^Tj) are expected to be obtained from the kinetic terms of the scalars, 
which are, however, not present in the effective Lagrangian obtained after a couple of manipulations. Therefore either 
that kinetic terms must be generated by quantum corrections or the correction terms to the Noether current should 
directly be calculated. What we have done in Sect. |l| is a calculation of the correction term directly concerned with 
the chiral symmetry breakdown. There we have calculated Fig. 2 and obtained the term proportional to dp^ip that we 
need now. That is, we have required that the diagram should satify 



SC 

Sdf^tp 



Sip = 2fd^ip{x) + . . 



(33) 
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which determines / in terms of 5](a;) after symmetry breaks down. Namely the Feynman diagram corresponding 
to Fig. 3 is now given by Fig.4, which gives massless particle (NG boson) contribution to the vertex. This is a 
multiplication of two diagrams, Fig. 2 and Fig.l and in between them there is a massless NG boson propagator which 
gives a pole term. 




Fig.4 

FIG. 4. Dynamical Nambu-Goldstone boson contribution to the axial vector vertex. 



Fig. 2 gives a factor — 2/g^, Fig.l does + (7/2) 75//, and the 1^ propagator does i/q^- After multiplying all these 
factors, we obtain — 2iE(p + g/2) 75 q^/q^ , which is exactly equal to the second term of the rhs of the WT identity 
Eq. (|0|). 

Finally we would like to describe another simple method to calculate the decay constant or /. In order to estimate 
Eq. ( p7| ) we lift global chiral symmetry to local one. That is, we introduce gauge fields associated with local chiral 
symmetry, calculate some diagram to obtain /, and then turn off local symmetry to have the original global chiral 
symmetry. The advantage to have local symmetry is to have a gauge coupling term with the Noether current, 
V!(a;) 7^75'(/'(a;), as 

4(x)V;(a:)7^75V'(a^), (34) 

which helps calculate /. We here neglect the chiral anomaly as well as its effects due to introduction of local chiral 
symmtery or A^^{x). This docs not affect the final results. The fcrmion-anti-fcrmion-NG boson vertex is given by 

f^(p + g,p) = j d^xd'ye-^'^^-^^y m{Al{xmy)m)\^).^^, (35) 

which corresponds to a bare vertex 17,^75. Here the subscript "amp" means vertex functions arc amputated. Taking 
an analogy to the ordinary elementary Higgs field, what we need is the following two-point function, 

V,{q) = I d'x e-'- (0| T (4(0)^(x)) |0),„p , (36) 

which is nothing but 2/ q^ in the limit of which will be shown below. With the help of Eq. (|3^) Eq. ( ^ ) is 

given by a fermion one-loop diagram depicted in Fig. 5 and equating this to 2/5^ gives the expression for / in terms 
of Yj{x). This is the simplest way to obtain the Pagels-Stokar formula compared with other methods. 




Pig.5 

FIG. 5. Fermion loop correction to the ylj^ — two point function. 



To see that the above two-point function Eq. (|3^) gives If q^, we expand the ordinary "gauge-invariant" Higgs 
kinetic term. 
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\D^{a'{x)+Mx))f 
^[{d^ ~ 2^A^'^) {a'{x) - z^(.t)) {8^ + 2*4) {<j'{x) + t^{x))] - 2/4^^ + 



(37) 



where the charge of th e scalar term a' {x)+iip{x) is -2 as shown by Eq. (32b) and use has been made of < 0|(T'(a:)|0 >= / 
as given by Eq. (|l6aD . Notice that a relative sign between 4V'7Ai75'0 a-nd 2fA^d^ip in Eq. ( ^T\ ) is plus which is the 
same as in Eq. (pll). 

There are some models in which there are still gauge fields remained after functionally integrating a couple of 
strongly coupled gauge fields. In this case there remains local symmetry and we have terms similar to Eq. ( ^ ) which 
give, as we have seen, masses to remained gauge fields, which will be a next problem to be solved. 
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